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Abstract 

The Nested Bethe Ansatz is generalized to open and independent boundary 
conditions depending on two continuous and two discrete free parameters. This 
is used to find the exact eigenvectors and eigenvalues of the A^-i vertex models 
and SU{n) spin chains with such boundary conditions. The solution is found for all 
diagonal families of solutions to the reflection equations in all possible combinations. 
The Bethe ansatz equations are used to find de first order finite size correction. 
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As is by now well known integrability in the bulk is a consequence of the Yang-Baxter 
equation (YBE) in two-dimensional lattice models and two-dimensional quantum field 
theory (QFT), (see for example ref. [0) . Recently interest is focushed in generalizing 
boundary conditions (b.c.) compatible with integrability. Periodic b. c. , twisted b. 
c. (by Cartan algebra operators) and K-matrix b. c. in Sklyanin's framework provide 
integrable models [0. Several new families of b. c. compatible with integrability has been 
recently found 0, ^, ^ . 

The quantum group invariant chains can be obtained by this method as special cases 
P, 1^, ^ 1^. In this construction the refiection equations appear as the new ingredient 
to the bulk Yang-Baxter equation when dealing with open boundary conditions. For the 
moment only models corresponding to rank one algebras has been solved with indepen- 



dent open boundary conditions 0, ^ and recently the t-J model |Tl[]. For algebras of 

(2) 

larger rank only the quantum group invariant case of the A2 , of the v4„_i and the t-J 



models has been solved ||9|, O, IT3 



In this letter we solve the integrable models found in 0. They are models with indepen- 
dent open boundary conditions associated to the algebras A„_i. We refer to the article 
V2\ for the notation and details of some parts of the solution (we will work with a spectral 



parameter shifted as^^^^^ + 7/2 with respect to [|T^ ) . 



The An-i algebras have non symmetric i?-matrices given in the hyperbolic regime by: 

iCM = 1 

1 < a,b < n. 

The diagonal solutions to the corresponding refiection equations were found in 0], 
they are: 



sinh ^_ 



'-/C(9) = 



sinh ^_ 



l<a<L 



+ 1 <a<n, (2) 



^ sinh(2^ + n7) sinh^+ - - ^ 
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^-KM = ^^"M2g + 7) e("-^"+^)^ sinh(e+ - g - n7) + i < , < (3) 

" sinh(2^ + 717) sinh^+ 

Notice that the quantum group invariant case is obtained in the hmit C,± 00. 
Using these solutions one can construct for an An-i algebra, a set of {n — 1)^ different 
transfer matrices depending on two arbitrary parameters by taking a '■+K~^ matrix and a 
'^~K~ matrix. We will call the transfer matrices constructed in this way '+'"t(^,u)), where 
u = {ujpo, are the first level inhomogeneities of a chain of length pq. 

The way to diagonalize these transfer matrices is parallel to the quantum group invariant 
case but we need som new ingredients. First we need to know the action of the doubled 
monodromy matrix operators on the reference state. The computation of this action is 
similar to the quantum group invariant case, but now we use the expression (in the 
quantum group invariant case K~{6) = 1). We have at level (k) of the nested Bethe 
ansatz construction: 



+ (A; _ 1)^/2) II 1 >= ^-K[''\e +{k- 1)7/2) II 1 >, 

+ 1)7/2) II 1>=0, 

+ ik- 1)7/2) II 1 >= A^\e + {k- 1)7/2)5,, II 1 >, 
B!t\9) ||1>^ 0, 

" ^ ' /=i sinh[^ + /if-') + fc7]sinh(e-/if-')+7) ' ^ 

^-Kr'\e + k^/2) = e-^7/2sinh[(e- + fc7/2)-g-A:7/2] k + l<d<l., 

sinh^_ ~ ~ 

+ A:7/2) = ,->^./2 ^^m- + kl/2) + e + k^/2] i_ + i < d < n, 

sinh^_ ~ ~ 

(4) 

where there are no summations over repeated indices and || 1 > is the reference state of 
the corresponding level. It is easily seen that the matrix ^- K^'^'^^\6 + /C7/2) obeys the 
reflection equation for a problem of weights R^^^^\0 + ^7/2) with indices from /c + 1 to 
n. Also the ^+K^{9) matrix obeys: 
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hKlt\e + ik- 1)7/2) = ^ik-i),/2^j^M^^J±fl [(„-fc+i)-2(d^fc+i)+ih 
^ ^ ' sinh(2^ + n7) 

sinh[(e+ - (fc - 1)7/2) + g + (fc - 1)7/2] A; < d < /+, 

sinh^+ ' — — ' 

^+K^^\9 + (k- 1)7/2) = c<^'^^^/^ smh(2g + fc7) ^[(^n-k+i)-2id-k+i)+ih 

sinh(26' + ^7) 

sinh[(^+ - (A; - 1)7/2) - g - (n - fc + 1)7 - (A; - 1)7/2] ^e+{n-k+ih+{k'iyi/2 ^ +i<d<n 

sinh^+ ' ~ ~ 

(5) 

The matrix ^+K^^\6 + {k — 1)7/2) obeys the reflection equation for a problem of weights 
R^^\e + {k- 1)7/2) with indices running from k to n. In the previous equation some 
unimportant factors have been kept to make this property clearer. 
Using eq.@, the transfer matrix can be written as: 



i+i 

d=l 



sinh[^+ + e-{n- 



sinh 



e 



Now we want to evaluate the action of the transfer matrix on the first level Bethe 
Ansatz wave function, given by: 



2<ij <n 

This is done using the commutation relations between the doubled monodromy matrix 
operators, which are the same as in the quantum group invariant case. The result is: 



'+'-t(0,cD)^ = 

sinh[e+ + g-(n-/+)7] sinh(e--g) ^ sinh(g + /xf ) sinh(g - /x^.^^ - 7) 

sinh ^+ sinh ^_ ^ f}^ smh{9 + ^if + 7) sinh(0 - ^if) 

sinh(2g) ^ sinh(g + u,) sinh(g - Uj) 
^ sinh(26' + 27) W sinh(6' + cu^ + 7) sinh(6' - cjj + 7) 
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^ sinh(g + + 27) sinh(g - /x}'^ + 7) 
fli sinh(^ + /ij.') + 7) smh(^ - iif ) 

+ u.t, 

where u.t are the unwanted terms. We are then led to a reduced problem of transfer matrix 
-with weights R^'^\9 + 7/2), reflection matrices given by ^+K^'^\6 + 7/2), 
^~K^'^\9 + 7/2) and inhomogeneities given by the first level Bethe Ansatz. The reduced 
eigenvalue problem is: 



^(1))X = '+'-A(2)(^; (9) 

Using that tr2 R^i2 ~^ 1 / '"'"-^2^'' (^ + 7/2) is a diagonal matrix is possible to see that 
the unwanted terms cancel if: 



z+«_^(2)/ _ sinh[e+ + /ij^^ - (n - /+)7] sinh(e- - /ij^^) 

' sinh^+sinh^_ 

-Q sinh(/xi^^ + /if ^) sinh(/ii^^ - /if ^ - 7) 



t/k sinh(/if ^ + /if ^ + 27) sinh(/if ^ - /if ^ + 7) 

sinh(/if ^ +Ui + -f) sinh(/if ^ - + 7) 1 / l / „ nn^ 
l\ sinh(/if)+..)sinh(/if)-..) ' 

The analytic properties of the eigenvalues can be also used as a short way to obtain 
these equations. One can follow this strategy level by level using equations When 
the Bethe ansatz level is mm(/+,/_) we can extract a common factor from the reduced 
problem of order mm(Z+,Z_) + 1. The same occur when the level is max(/+,/_). In 
this case the reduced problem has the matrices of the quantum group invariant case. 
These factors are: 



fi-, k = L 



sinh(^_ + 6* + /_7) 

sinh(^__) 
sinh(^+ -9-n-f) 

sinh(^+) 

sinh(^_i_ — 9 — nj) sinh(^_ + 9 + I'j) 
sinh(^+) sinh(^_) 



fh k = l^ 
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where k is the Bethe ansatz leveL Using the previous resuhs the recurrence relation for 
the eigenvalues can be obtained. The general expression is: 



^(,) |i sinh[g + + {k- 1)7] sinh(g - /if ^ - 7) 

/ii sinh(^ + /if) + k^) smh{e - /if)) 

sinh[2g+ (fc - 1)7] smh[e + /i^^"^) + (fc - 1)7) sinh(g - /i^^-^)) 
^ ^ ^sinh[2^ + + 1)7] jj^ smh{e + /i(*^-i) + kj) smh{e - /iC^-i) + 7) 



^ sinh[^ + /if) + (A; + 1)7] sinh(^ - /if) + 7) 
jii sinh(^ + /if) + kj) smh{9 - /if)) 



1 < A; < n - 1, /if) = cj,-, A(")(^, /i("-i)) = 1. (12) 
The expressions of the functions a^^\6)^h^'^\6) are given by: 

mm(/_|_, /_) — 1, 

a('<)(^)=a, b^^<\e) = fi^, k = mm{l+,l.) ■.= U, 

a^^\e)=gi^, 6W(^) = 1, /< + l < A;<max(/+,/„)-l :=/>-l, 

a^'>\6)=g,^, b^^>){e) = fi^, k = l^, 

aW(^) = 6W(^) = 1, /> + l<A;<n-l, (13) 
where in the case 1+ = l_ = I, fi^ = fi^ = fi and gi^ = a. In the previous expressions: 



sinh[^4. + 9 — {n — l+)'y] sinh(^_ — 
sinh ^+ sinh ^_ 



a = ^ ^ ^ ^ e("-'+)^, 



^ sinh[e+ + e-{n- 
^ sinh ^+ ' 

smh 4- 



We can now use the recurrence formula (|12|), and with the help of equations ([T3|,[T 
find the expression for the eigenvalue of '+'~t(^, cD). The result is: 



A//1 ~N -TT sinh(6' + Co",) sinh(6' - u;,- 



j-J-^ sinh(6' + iUj + 7) sinh(6' — luj + 7) 
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^ (fc) smh(2g) sinh(2g + 7) 



n 



smh[^ + /if + k-f] smh{e - /if "'^ + 7) 



1 sinh[0 + fif~'> + {k- 1)7] sinh(0 - fi 



smh[9 + /if ^ + {k- 1)7] sinh(^ - /if ^ - 7) 

] si: 

where the product over p„ is equal to one and /if ^ = ooj. The value of g'^^\9) is given by 



/ii sinh[e + /if + A;7]sinh(^-/if ) 



9^''\0) = fiji_, ly + l<k<n. (16) 
The expression for the eigenvalue ([l5| , p!6D gives the quantum group invariant result in 



the limit ^± ^ 00 |[T^. For the SU{2) case this formulas reduce to those in reference 
0. One can see that the "mixed" cases where /+ 7^ /_ give imaginary eigenvalues in the 
trigonometric regime. For the cases when = l_ = I the eigenvalues are essentially real, 
we can get rid of the exponential factor by a redefinition of the the reflection matrices 
— > e~^^"^''^^ 1 11]. This factor has been kept to make an easier connection with the 



quantum group invariant case. 

The Bethe ansatz equations that the roots /if ^ have to obey are: 

^(,). (,) sinh[/if ) + /if) + {k- 1)7] sinh(/if ) - /if) - 7) _ 
sinh[/if') + /if) + ik + 1)7] sinh(/if ) - fif + 7) 

sinh(/if ) + /if + fc7) sinh(/if ) - /if - 7) 
M sinh[/if ) + /if + (fc + 1)7] sinh(/if ) - /i('+')^ 

n 



^''-^ sinh[/if ) + fif-'^ + {k- 1)7] sinh(/if ) - /if '^)) 
/Ji sinh(/if ) + /if ~"^) + k^) sinh(/if ) - /if ""^^ + 7) 

l<A;<n-l, l<i<Pk. (17) 



In the previous expression the function h^'^^ [6) is given by: 

smh(e_+/ir) + Z_7) 
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h = l- = l. h^Hf^?) 



sinh(,^+ — /if"'"'' — 727) 

1 



From this last expressions we see that the eigenvalues of the transfer matrix in the 
trigonometric regime when = = / are essentially real and otherwise imaginary. 

Let us now derive the solution of the present NBA equations in the thermodynamic 
limit for the ground state and in the gapless regime. For this purpose, it is convenient to 
relate these equations with those of the periodic case by means of the change of variables, 
(see [161 ra): 



1 < S < Pfc 



-v. 



(k) 



1< k < n - 1, 



where i/]'^^ 



(k) 



ik'~f/2. To pass to the gapless regime we make also the changes 



(k) 



IV. 



(fc) 



. Then, following standard techniques one introduces 



a density of roots at every NBA level (from now on po = ^)'- 



lim 



As explained in refs. P, a hole at A = shows up in the ground state. Then 



is a regular and continuous function. In the N ^ 00 limit, eqs.(p!7|) become through 
standard techniques a system of integral equations for the densities cr'(A). 



a\X)-Y. / rf/ii^fcm(A-/x)a'"(/x) 
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^ $'(A,7/2)-;^$'(A,(vr-7)/2) 



That is, 



+^ <|.'(A, 7/2) - - E + ^ xUA) (19) 

Xfc(A) = -^log/iW(A) (20) 



(A) = 2tX + <|.(A, + |/-) , Xi^ (A) = -2^A - $(A, e+ - 7(r^ - y )) 
Xfc = when A; ^ /+, /_. (21) 
When /+ = /_ = /, x/(A) =Xz-(A) + xz+(A) (22) 



Here 



$(^,7) = i log 



sinh(27 + z) 
sinh(i7 — z) 



When Xfe(A) = we recover the SUq{n) invariant case solved in ref. [12|. The density 
of real roots in our case results equal to cr'(A) in ref. [T^l plus the term 



5a\k) 



-oo 2ti 



5a\k) 



(23) 



1 



sinh/c(- - r]^)Rii_{k) - sinh/c(- - r]+)Ra_^{k) 



sinh fc| 



2i 



[i?,.(0)-i?„40) 

One can see that in the quantum group invariant case this term goes to zero. Here Riv{k) 
stands for the resolvent kernel: 



Rii'{2x) 



sinh(7ra;) sinh[7x(?T, — />)] sinh(7x/< 
sinh[x(7r — 7)] sinh(7x?T,) sinh(7x) 



and 



^ 2 



(We assume vr > ri± > 7/2). We will see what is the change in the free energy for the 
case when /+ = /_= Z, in the other cases one obtains imaginary contributions due to the 
fact that the eigenvalues are imaginary, we also perform the redefinition of the matrix 
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mentioned above. The change on the free energy compared to the SU{n)q invariant case 
is: 



cj-rn\ 4 f°° dk sinhfc6' sinh[A;7(ra — /)/2] sinh[A;(r7+ — ?7_)/2] cosh[A;(7r — 77+ — ?7_)/2] 
df{9) - 




N Jo k sinh (7^/0/2) sinh (7rA;/2) 

sin[e+ + e-{n- /)7] \ 1 f sin(e- - 9) \ 
sin^+ J ^ 1 sin ^_ J 

This term represents the surface energy induced by our general boundary conditions 
as a function of ^+ and The derivative of Sf{6) at ^ = times — provides the 
surface change on the ground state of the associated magnetic chain. 
In conclusion we have obtained the expression for the eigenvalues and nested Bethe ansatz 
equations for the open A„_i models of reference 0. For some regimes the "mixed" cases 
give complex eigenvalues. For the cases where /+ = /_ = / the eigenvalues are real 
for all the regimes even though the hamiltonians are not hermitian. We have obtained 
the first order finite size corrections in this case. It would be interesting to study the 
thermodynamic properties of these chains. In the other hand these hamiltonians, even in 
the mixed cases, could have some application in the framework developed in for the 
reaction diffusion equations. 



AGR would like to thank L. A. Ibort for comments. 
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